
cd

use "Main.dta", clear

* density and summary stats for elasticity estimates

twoway kdensity ne20, lcolor(black) lpattern(solid) || kdensity le20, lcolor(black) lpattern(dash) || kdensity se20, lcolor(black) lpattern(shortdash) legend(pos(1) ring(0) col(1) lab(1 "NE20") lab(2 "LE20") lab(3 "SE20"))
tabstat ne20 le20 se20, stat(mean med sd) by(city) col(stat) long

* correlations between ne20 and census controls

corr ne20 log_area log_pop urban nwnp age1844 school1620 mfgest_pc mfgwork log_mfgva farms_pc acres log_avgcrop log_avgvalue homeownclr
corr ne20 log_area log_pop urban nwnp age1844 school1620 mfgest_pc mfgwork log_mfgva farms_pc acres log_avgcrop log_avgvalue homeownclr if city==1
corr ne20 log_area log_pop urban nwnp age1844 school1620 mfgest_pc mfgwork log_mfgva farms_pc acres log_avgcrop log_avgvalue homeownclr if city==0
